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Abstract: Different methods for retrieving atmospheric profiles in the presence of clouds from 
hyperspectral satellite remote sensing data will be described. We will present results from the 
JPSS cloud-clearing algorithm and NASA Langley cloud retrieval algorithm. 

OCIS codes: [010.1280] Atmospheric composition; [010.5620] Radiative transfer 

1. Introduction 

Hyperspectral satellite sensors such as Infrared Atmospheric Sounding Interferometer (IASI) and Cross-track 
Infrared Sounder (CrIS) have high spectral and spatial resolution. They provide abundant information on 
atmospheric and surface properties. In order to retrieve atmospheric temperature, water, and trace gas vertical 
profiles from these high spectral resolution data, one has to account for cloud contributions to the top of atmospheric 
(TO A) radiance spectra. Usually there are two methods for dealing with clouds: cloud-clearing (CC) and cloud- 
retrieval (CR)[1,2,3,4]. The CC is used by the AIRS level 2 data processing algorithm and by the Joint Polar 
Satellite System (JPSS) Cross-track Infrared and Microwave Sounder Suite (CrIMSS) algorithm. At NASA Langley 
Research Center (LaRC) [5, 6, 7, 8, 9], we have developed a retrieval algorithm, which explicitly retrieves cloud 
properties together with other parameters such as atmospheric temperature, moisture, and trace gases profiles, 
surface skin temperature and emissivity. We will present results of testing the CrIMSS Environmental Data Record 
(EDR) operational algorithm and the LaRC CR method using IASI satellite data. 

2. JPSS CrIMSS algorithm 

The operational CrIMSS EDR an algorithm was originally designed to run on large IBM computers with dedicated 
data management subsystem (DMS). We have ported the operational code to simple Linux systems by replacing 
DMS with appropriate interfaces. The ported code is capable of reading data from both the operational code and 
from the CrIMSS science code. The detail of the CrIMSS EDR algorithm is described in reference [1]. The cloud- 
cleared radiance ^ dr for channel i can be expressed as a linear combination of the measured radiances 
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Where rj x ...rj K are unknown channel-independent constants, % % are the measured radiances and at least K+l 

field of views are needed to solve for K cloud formations. We will test the CRIMSS algorithm performance by 
generating proxy data from the IASI observations. The details will be given in the result section. 

3. NASA Langley cloud retrieval algorithm 

One of the uniqueness of the LaRC CR algorithm is that it uses a principal component based radiative transfer 
(PCRTM), which is capable of using all IASI spectral channels for the retrieval. It predicts projection coefficients of 
channel radiance onto a set of orthogonal eigenvectors (PCs), which capture the spectrum information in a 
descending order according to the associated eigenvalues. The channel radiance can be easily generated by a linear 
combination of PCs. The linear coefficients are called PC scores (or projection coefficients). Generally, only 100- 
200 PCs are enough to regenerate the channel radiance to required accuracy. Since the dimension of the 
hyperspectral data is reduced from thousands to ~100, the computational cost for both forward radiative transfer 
calculations and inversions can be reduced dramatically. The inversion algorithm is based on a non-linear 
Levenberg-Marquardt method with climatology covariance and a priori information as constraints: 
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Where x represents state vector, the subscripts n and a priori, respectively. is the Jacobian matrix. Y m is the 
transformed observations (in this case an EOF transformed radiance spectrum). Y n is the forward model calculated 
PC scores using the state vector obtained from the nth iteration. S y and S a are error covariance matrices associated 
with observation (including forward model errors) and background state vector, respectively. All relevant parametes 
such as cloud properties, atmospheric profiles, and surface properties are retrieved simultaneously. The spatial 
resolution of the LaRC approach is higher than the CrIMSS EDR algorithm since it does not need to use multiple 
FOVs to perform cloud clearing. 

4. Results 

The IASI on board the METOP-A satellite provides all the information needed to generate CrIS proxy data. The 
IASI instrument is a Fourier Transform Spectrometer (FTS) with a spectral sample interval of 0.25 cm" 1 and a 
continuous spectral coverage from 645 to 2760 cm' 1 . The CrIS instrument has three spectral bands with spectral 
coverage of 648.75 - 1096.25 cm" 1 , 1207.5 - 1752.5 cm" 1 , and 2150.0 - 2555.0 cm" 1 , respectively [5]. The 
transformation from an IASI radiance spectrum to a CrIS proxy spectrum is a rigorous mathematical operation. The 
generated proxy data have been successfully ingested into the CrIMSS OPS code. This data set is very useful for 
algorithm testing and validation because it provides realistic scenes (atmosphere, surface and cloud), which are 
based on observations. It will help us identify areas of deficiencies in both forward and inverse models, testing 
robustness of the operational code, check error handling capability of the OPS code. It will be used for bias 
characterization and algorithm tuning, and EDR performance evaluations. Details of the results of this study will be 
presented at the conference. 

We have applied the PCRTM radiative transfer model and inversion algorithm to hyperspectral sensors such as 
AIRS, IASI and NAST-I. The errors between the PCRTM calculated brightness temperature spectra and those 
calculated using line-by-line radiative transfer model are typically less than 0.05 K, a value much smaller than the 
instrument noise. Figure 1 shows two IASI brightness temperature spectra, one observed by the IASI (red) and the 
other calculated using the PCRTM (blue). The cloud properties are retrieved using the method described in the 
previous section. The retrieval algorithm identifies the cloud phase as ice. The retrieved cloud top pressure is at 273 
hPa. The cloud effective size is 38 mm and the cloud visible optical depth is 0.462. The bottom panel of the figure 
shows the difference between the observed and the calculated IASI spectra. They agree with each other within the 
IASI instrument noise levels (shown in red lines) in most of the spectral regions. 
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Figure 1. Top: Observed IASI spectrum and the PCRTM calculated spectrum. Bottom panel: Difference between observed and 
calculated IASI spectra (blue curve) and the IASI instrument noise converted to brightness temperature unit (red curves). 

Figure 2 shows retrieved global monthly mean (July 2009) temperature at 300 mbar (top left), surface skin 
temperature (top right), surface emissivity at 1140 cm" 1 (bottom left), and carbon monoxide mixing ratio near 500 
mbar. The global distributions of atmospheric temperature and surface skin temperature agree well the ECMWF 
monthly means. The temperatures are high in the tropical region and low in polar region with regional and seasonal 


variations well captured by the retrieval system. Studies done in the past have shown that the retrieved temperature 
and moisture profiles agree well with collocated radiosondes. The surface emissivity retrieval provides emissivity 
spectrum at each IASI footprint when the cloud is not thick. The map shown in this figure clearly displays the low 
emissivity values of sand features in desert regions. The CO mixing ratios are elevated near either pollution sources 
(southeast Asia) or biomass burning sources (central Africa). 



Figure 2. Retrieved global monthly mean (July 2009) temperature at 300 mbar (top left), surface skin temperature (top right), 
surface emissivity at 1 140 cm-1 (bottom left), and carbon monoxide mixing ratio near 500 mbar. 

4. Conclusions 

In preparation for the launch of NPP satellite, we have tested the CrIMSS operational code using proxy data 
generated from the IASI spectra. The CrIMSS EDR algorithm works well with the proxy data. We have also tested 
another algorithm developed at NASA Langley Research Center to retrieve atmospheric profiles along with cloud 
parameters. The method has been successfully applied to the IASI spectra directly. 
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